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Random tensor models for a generic complex tensor generalize matrix models in arbitrary dimensions 

and yield a theory of random geometries. They support a 1/N expansion dominated by graphs of 

spherical topology. Their Schwinger Dyson equations, generalizing the loop equations of matrix 

fvj . models, translate into constraints satisfied by the partition function. The constraints have been 

shown, in the large N limit, to close a Lie algebra indexed by colored rooted _D-ary trees yielding 

("—^ ' a first generalization of the Virasoro algebra in arbitrary dimensions. In this paper we complete 

rvq , the Schwinger Dyson equations and the associated algebra at all orders in 1/A^. The full algebra 

of constraints is indexed by D-colored graphs, and the leading order D-ary tree algebra is a Lie 

j^ ' subalgebra of the full constraints algebra. 
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I. INTRODUCTION 



O Random matrix models encode a theory of random Riemann surfaces dual to the ribbon Feynman graphs generated 

D ■ by their perturbative expansion [l|. The amplitudes of the Feynman graphs depend on the size of the matrix, TV, 

i-G I and the perturbative series can be reorganized in powers of 1/N [4]. At leading order only planar graphs contribute 

[3| . The planar family is summable and undergoes a phase transition to continuous Riemann surfaces [^, |5| when 

the coupling constant approaches some critical value. Single and multi-matrix models have been very successful in 

^ ■ describing the critical behavior of two-dimensional statistical models on random geometries [6|-[9| and via the KPZ 

•/"J ' correspondence nQ-[ lt3 on fixed geometries. Sublcading terms in the 1/iV expansion can be access through double 

^ ! scaling limits J14|-|l6t. 

^1 ' Random matrices are analyzed with many techniques, like gauge fixing to the eigenvalues of the matrix, orthogonal 

. ] polynomials and so on [l|. Of particular interest for this paper are the Schwinger Dyson equations (SDE) (or loop 
C^ . equation) |17h2(J | which translate into a set of Virasoro constraints satisfied by the partition function. 
^— ^ I The success of matrix models in describing random two dimensional surfaces inspired their generalization in higher 

~_" , dimensions to random tensor models [2l| - (23| (see also |24| - [30| for more recent developments). The corresponding 
theory of random higher dimensional geometries was initially hoped to give insights into conformal field theory, 
statistical models in random geometry and quantum gravity in three and four dimensions. In spite of these initial 
high hopes, tensor models have for a long time been unsuccessful in providing an analytically controlled theory of 
random geometries: until recently all the nice aspects of matrix models could not be generalized to tensors, as their 
d ■ 1/iV expansion was missing. 

The situation has changed with the discovery of colored |3ll-l33| rank D > 3 random tensor modeltQ. Their 
perturbation series supports a 1/iV expansion |35l437l |. indexed by the degree, a positive integer which plays in higher 
dimensions the role of the genus. Wc emphasize that the degree is not a topological invariant. Leading order graphs, 
baptized melonic [38[ , triangulate the D-dimensional sphere in any dimension |35l. |36| . They form a summable family 
as they map to colored rooted D-ary trees [3^. Colored random tensors [33 gave the first analytically accessible 
theory of random geometries in three and more dimensions and became a rapidly expanding field [iol - liq , (see also 
[47| for some related developments). In particular the first applications of random tensors to statistical models in 
random geometries have been explored [48l - l5]| . 

The results obtained for the colored models have been shown to hold in fact for tensor models for a single, complex, 
generic (that is non symmetric) tensor [52l | . The colors become a canonical bookkeeping device tracking the indices of 
the tensor. The theory thus obtained is universal |53J and constitutes the only analytically controlled generalization 
of matrix models to higher dimensions wc have so far. 
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^ In D = 2 colored matrix models have been studied [34|. but the colors do not play the key role they play in three and more dimensions. 



A posteriori one understands why the discovery of the \/N expansion for tensor models has taken so fong: in order 
to have any control on the amplitude of graphs one needs to be able to track the indices of the tensors (this is achieved 
by the colors). This is impossible if one uses tensors with symmetry properties under permutations of the indices. 
Tensor models for generic tensors have been generalized to tensor field theories, |54| - [56| in which the quadratic part 
is not invariant under the ®^U{N) symmetry. The non trivial quadratic part generates a renormalization group flow 
and the 1/N expansion is recovered dynamicall y. W e already possess an example of a tensor field theory which is 
asymptotically free |55| . Due to the universality [53| of tensor measures one is tempted to conjecture that asymptotic 
freedom will be a feature of all tensor field theories. 

In order to study this universal theory of random tensors one must generalize to higher dimensions as many of 
the tools which were so useful in matrix models as possible. This turns out again to be a non trivial problem. 
For instance, although various generalizations of eigenvalues to tensors have been proposed, they lack most of the 
interesting properties of the eigenvalues of matrices (for one the tensor does not diagonalize under a change of basis). 
A generalization of the notion of determinant to tensors (formats), the Gel'fand hyperdeterminant [57| exist (and 
interestingly enough it requires that the tensor indices are distinguished), however one does not possess an explicit 
formula for it (except for the (1,1,1) format, in which case one gets just the well known Cayley hyperdeterminant). 
However, 2 is a rather small value for a supposedly large N. 

The only set of techniques which can be generalized relatively straightforwardly to tensors are the SDE's. At leading 
order the SDE's of tensor models have been derived in [53|. As one writes a SDE for every observable of the model, 
the equations are indexed by the observables. At leading order only melonic observables contribute and, as they 
map to trees, the SDE's and associated Lie algebra of constraints are indexed by trees. The definition of the algebra 
relies on a gluing operation for trees, 7i •v 72, which generalizes the addition of integers (in terms of composition of 
observables for matrix models the addition encodes the gluing of two cycles of lengths p and q) to an "addition" of 
observables indexed by trees. Thus the Virasoro (strictu sensu the positive part of the Witt) algebra of constrains of 
matrix models becomes 

[Lm,L„] = {in-n)L„i+n =^ \C-r^,llr^ = X! ^71*^X2 - X! ^Tl^vTx ■ (1-1) 

In [5^ it is proved that the constrains hold at leading order in N ^ that is limjv->oo A^ '^^I^'Ct-Z = 0, with Z the 
partition function of tensor models. In this paper we complete the algebra of constrains to all order in 1/A. This 
time the observables (hence the SDE's and the algebra of constraints) arc indexed by ZJ-colored graphs B with a 
distinguished vertex v. The observables can be glued by a ^^(^^r,') operation and we obtain a Lie algebra of constraints 
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obeyed at all orders in N . When restricted to melonic observables the *(i,.5) composition of observables particularizes 
to the -ky composition of the associated trees. The melonic observables are closed under this gluing. We thus not only 
dispose of the full set of constraints for random tensor models, but, non trivially, the D-ary tree subalgebra closes 
into a Lie subalgebra of the full constraints algebra at all orders in 1/N. 

The algebra of constraints is the starting point of the study of the continuum limit [52l of tensor models. In matrix 
models the continuum theory, Liouville gravity coupled with various matter fields at a conformal point, is identified 
using the theory of (unitary and non unitary) representations of the Virasoro algebra. The continuum operators, 
which are composite operators in terms of the loop observables, are hard to identify Q. One uses the theory of 
representations of the Virasoro algebra to infer the appropriate correspondence. The same must be done in higher 
dimensions. We need first to study and classify the central extensions and unitary representations of this algebra. 
Once the field content of the continuum theory is identified, we need to translate the continuum observables in terms 
of the observables of tensor models. This study will provide a landscape of universality classes of continuum theories 
of random geometries accessible using tensor models in arbitrary dimensions. The algebra of constraints identified in 
[58| and completed in this paper at all order constitutes the first step towards applying tensor models to the study of 
conformal field theories and quantum gravity in arbitrary dimensions. 

This paper is organized as follows. In section |ll] we recall the generic one tensor models and their 1/N expansion. 
In section [TTTl we define two graphical operations on arbitrary observables of such models, derive the SDE's and prove 
that they close a Lie algebra. 



II. THE 1/iV EXPANSION OF TENSOR MODELS 
A. Tensor invariants and action 



• Hd be complex vector spaces of dimensions Ni , ■ 



■ Nd . A covariant tensor T„ 



of rank D (or 



Let Hi, 
a (A''i — 1, . . . No — 1) format) is a collection of Hili ^i complex numbers supplemented with the requirement of 
covariance under base change. We consider tensors T transforming under the external tensor product of fundamental 
representations of the unitary group ®f^-JJ{Ni), that is each U{Ni) acts independently on its corresponding Hi. The 
complex conjugate tensor T„-^...nD is then a rank D contravariant tensor. They transform as 
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(2.3) 



We will denote the indices of the complex conjugated tensor with a bar, and use the shorthand notation n = 
(ni, . . . , Ud)- We restrict from now on to Ni = N, Mi. 

Among the invariants built out of T and T we will deal in the sequel exclusively with trace invariants. They are 
obtained by contracting two by two covariant with contravariant indices in a polynomial in the tensor entries, 



Tr(T,r) = ^]J(5„i,fii Tn,......Tfi^... 



(2.4) 



where all indices are saturated. A trace invariant has the same number of T and T. As T„j^...„j-, transforms as 
a complex vector under the action of the unitary group on only one index, one can use the fundamental theorem 
of classical invariants for U{N) (whose origins can be traced back to Gordan [531, see [60| and references therein) 
successively for each index and conclude that all invariant polynomials in the tensor entries write as linear combination 
of trace invariants. 

Trace invariants are labeled by graphs. To draw the graph associated to a trace invariant we represent every T by 
a white vertex v and every T by a black vertex v. We promote the position of an index to a color: ni has color 1, 
712 has color 2 and so on. The contraction of two indices Ui and fii of two tensors is represented by a line P = {v, v) 
connecting the corresponding two vertices. Lines inherit the color of the index, and always connect a black and a 
white vertex. Any trace invariant is thus represented by a closed D-colored graph B [39[. The graph B is bipartite 
(it has black and white vertices denoted v and w, such that all lines connect a black and a white vertex) and its lines 
P = (ii, v) have a color i = 1, 2, . . . D. All the vertices of B have coordination D with all the edges incident to a given 
vertex having distinct colors. Some trace invariants for rank 3 tensors are represented in figure [TJ 





Figure 1. Graphical representation of trace invariants. 



The trace invariant associated to the graph B writes as 



D 



TMT,f)= Yl ^{n^M-n U^n^T^^^, with -5f.„_..j=n U S„.,. 



(2.5) 



where P runs over all the lines of color i of B. The product of Kronecker delta's encoding the index contractions of 
the observable associated to the graph B, S?^^, ^j,, is called the trace invariant operator with associated graph B J5q . 

Trace invariant operators factor over the connected components of the graph. We denote r2j. the set of _D-colorcd, 
connected graphs with 2fc unlabeled vertices. 

Formally, a D colored graph with labeled vertices is defined by an incidence matrix e^,g whose entries are {ii, . . -ik} 
if the vertices v and v are connected by lines of colors ii, . . . , z^. In particular e^y = if w and v are not connected by 
any line. An element B g Fj^. is an equivalence class of incidence matrices related to one another by permutations 
of lines and columns (corresponding to relabellings of the vertices) . One can write the trace invariant associated to B 
directly in terms of the incidence matrix of any representative graph with labeled vertices Yiv v Tii^e - '^""n"- 



The subgraphs with two colors of a D-colorcd graph B are called faces. We denote the number of faces of a graph 
B by -Fg. They will play an important role in the next section. The graphs with D = 3 colors represented in figure [T] 
posses three type of faces, given by the subgraphs with lines of colors (1, 2), ( 1,3) and (2, 3). 

To every graph B we can associate a non negative integer, its degree uiiB) |36l. 1371. |39| . The main property of the 
degree is that it provides a counting of the number of faces Fq: for a closed, connected graph with D colors and 2pB 
vertices the total number of faces is [33, |3^, [52| 

Fs = (^^-lK£^p^ + {D-l)- ^^^ o^iB) . (2.6) 

The degree provides in higher dimensions a generalization of the genus of ribbon graphs, and indexes their the 1/A'^ 
expansion. It is however not a topological invariant, but it combines topological and combinatorial information about 
the graph [4^. Of course the degree can be defined for graphs with D + 1 colors (say 0, 1, ... I?) , and the number of 

faces of graph with D + 1 colors is Fg ~ ' ^ ' pg + D — {d^i)\ ^(S)- Another important property of the degree 
(see for instance [33]) is the following. Consider a closed connected graph Q with D + 1 colors 0, 1, ... I?, and denote 
Bp its connected subgraphs of colors 1, . . .D (where p = 1 . . .\p\ labels the connected components). Then 

M 

io{g) > Dj2^iBp) ■ (2.7) 

p=i 

Going back to invariants one can build out of a complex tensor, we note that there exists a unique D-colored graph 
with two vertices, namely the graph in which all the lines connect the two vertices. We call it the ZJ-dipole and its 
associated invariant is 

D 



Trdipoic(T, f)=Y,Tn Ts [n ^n,",] ■ (2-8) 
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The most general single trace invariant action for a non-symmetric tensor is 

oo 

5(r,T)=iiTrdipoic(T,r) + ;^ Y, tBTi-B{T,f), (2.9) 

where tg are the coupling constants associated to the _D-colored graphs B and we singled out the quadratic part 
corresponding to the I?-dipole. In equation ()2.9p one sometimes adds a scaling factor iV~<°-2''"' ' for every trace 
invariant (as a;(B) > this suppresses some of them). Adding this extra scaling would simplify some formulae but 
does not modify anything in the sequel. We will treat in this paper the most general single trace rank D tensor model 
defined by the partition function 

Z(te) = cxp{-F{tB)) = f dfdT exp (-N^-^S{T,T)) . (2.10) 

B. Graph amplitudes 

The invariant observables are the trace invariants represented by D-colored graphs. The Feynman graphs con- 
tributing to the expectation of an observable are obtained by Taylor expanding with respect to tg and evaluating 
the Gaussian integral in terms of Wick contractions. The D-colored graphs B associated to the invariants Tiig{T,T) 
in the action act as effective Feynman vertices. The effective vertices are connected by effective propagators (Wick 
contractions, pairings of T's and T's). A Wick contraction of two tensor entries Ta-^,,,ao 3'iid Tpj.pj^ with the quadratic 
part (|2.8p consists in replacing them by j^d^i^ Y[i=i ^^iPi ■ ^'^ '^ill represent the Wick contractions as dashed lines 
labeled by the fictitious color 0. The dashed lines of color are thus very different from the solid lines of colors 
1, 2, . . . _D: they identify all the indices of the two vertices (one white corresponding to T and one black corresponding 
to T) it connects (recall that the lines of colors 1, 2, ... D identify only one index each). An example of a Feynman 
graph is presented in figured) 

The Feynman graphs Q are therefore {D + l)-colored graphs. We will keep the notation B for the D-colored graphs 
and denote G the {D + l)-colored graph. We call the connected components with colors 1, . . .D oi Q, denoted Bp with 
p = 1, . . .\p\, the D-bubbles of Q. For instance the 3-bubbles with colors 123 of the graph in figuredjare the subgraphs 
made of solid lines. 
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Figure 2. A Feynman graph of tensor models. 



A graph Q has two kinds of faces: those with colors i,j = !,...,£) (which belong also to some D-bubble i3), and 
those with colors 0, i, for i = 1, . . . , D. We denoted F'g the number of faces of colors Oi of Q. 
The free energy has an expansion in closed, connected (D + l)-colored graphs, 

where s(t/) is a symmetry factor and \p\ is the number of effective vertices (that is subgraphs with colors 1 ... D or 
D-bubbles). The amplitudes of Q is 

p {«",««} p ;°=(-u,o)ee « 

An index n^ is identified along the lines of color i in Bp and along the dashed lines of color 0. We thus obtain a free 
sum per face of colors Oi, so that 

A{Q) = Hp!^ iv(c-i)lp|-P-i)l''VE, 14"! , (2.13) 

xr I 

where |/°| denotes the total number of lines of color of Q. But ^^ Fg^ ~ Fg — ^ i^p, where Fg denotes the total 
number of faces of Q and Fp the number of faces of the D-bubble Bp. Using (|2.6p for each Bp and for Q (taking into 
account that Q has 13 + 1 colors) and noting that \l^\ = p, with p the half number of vertices of Q, we obtain 
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A(6;) = ^^p^'^' ^-D-jjj^^(e)+^,j^E|,'ii'<^(Sp) 



fi^ TV "' ' ' "'■"-^'■■K '- -^=^ ' "7 , (2.14) 

with Li!{Q) the degree of the graph Q and ci;(;Bp) the degree of Bp. The amplitude of a graph tj is thus at most A^^ and 
it is suppressed with the degree uj{G). This is the 1/iV expansion for random tensor models [33, 112, IS^l ■ Expectation 
values of the observables have similar expansions in 1/N 

TrB,(r,r)...TrB,„,(r,r)^^ ^ ^c_i„i(B_i)__i_^(c;)+_^5^iH^^(g^)rip£a^ ^ ^2.15) 

eDu„B„ 1 

where the sum runs over all _D + 1 colored graphs Q with |a| marked D -bubbles Ba, and Bp denotes all the D-bubbles 
of Q (hence the a's arc some of the p's). 

III. SCHWINGER DYSON EQUATIONS 

The Schwinger Dyson equations of the model write in terms of two graphical operations. They encode the changes 
an observable B undergoes when adding a line of color 0. The /^-colored graphs B represent D ~ 1 dimensional closed 
connected pscudomanifolds [Sj] . When seen as subgraphs of some D + 1 colored graph Q they become chunks of a 



D dimensional space. The chunk associated to B is the topological cone over the pseudomanifold [S^l- The same 
holds for matrix models: the invariants Tr[(MM '^)''] are either seen as loop observables, or as polygons when acting 
as effective Feynman vertices. The polygon represented by Tr[(MM''')''] when it appears as an effective vertex in a 
Feynman graph is the cone over the loop. 

The usual loop equations of matrix models encode the change of the loop observables by the addition of a line: a 
loop can either merge with another loop or it can split into two loops. In the dual triangulation this is seen as the 
gluing of two segments bounding loops. Either the two segments belong to two distinct loops, and then one obtains a 
larger loop, or they belong to the same loop, in which case the loop splits into two. A loop observable Tr[(MM'f)'?+^] 
is fully characterized by its "length" q. The merging of two loops of length p and q leads to a loop with length p + q: 
the merging is just a graphical representation of the addition of integers. The SDE's of the model are indexed by 
loops, hence by integers. They can be transformed into a set of constraints on the partition function LpZ = 0. The 
latter respect the Virasoro algebra (l7l - l2q |. [Lp, L^ = {p — (l)Lp+q. Graphically this equation is understood as follows. 
Given two loops of length p and g, in how many ways can I compose them? I can either glue the loop of length q on 
the loop of length p (and one has p+1 choices for where to insert the second loop) , or I can glue the loop of length p 
on the loop of length q (and one has <? + 1 choices for where to perform this insertion) In all cases I obtain a loop of 
length p + q, associated to the operator Lp+g. 

A similar picture holds, with the appropriate generalizations, in higher dimensions. 

A. Graph operations 

The colored gluing of two graphs. Consider two D colored graphs Si and B2 and chose two vertices vi £ Bi and 
V2 G B2. We define the colored gluing of graphs at vi and V2 as the graph Bi *(i,i,c2) ^2 obtained by deleting the 
vertices vi and V2 and joining all the lines touching them pairwise respecting the colorings. 

This operation can be performed in two steps. Consider the trace invariants associated to Bi and B2 and connect 
the two vertices vi and V2 by a dashed line of color 0. Forgetting for a second the scalings with N , this line will 
identify all the indices of the tensor associated to vi with the ones of the tensor associated to V2 

E ^s^.n(^^np<0^g^^n= E (n n M(n^«p<0(n n ^<^j 

= ,5fr<;if '"'^ , (3.16) 

as vi is a vertex in Bi and V2 is a vertex in B2- The gluing is represented graphically in figure [3l 




'V^, 




Figure 3. Graphical representation of the colored gluing of two graphs. 

Of course the gluing preserves the colorability, hence Bi -^(vi,v2) ^2 is a connected D-colored graph with colors 
1, . . . D. At the level of the incidence matrix, one builds from e^^ and e„| the incidence matrix 

(e^.l,iiv,veBi 
v^vi,v^V2 e^,^'''^'-'''^' = L^',iiv,veB2 =e^^Ue^?u[e^^^ne%). (3.17) 

[e^^,ne% xiveB2,veBi 

The colored contraction of a graph. The second operation is similar to the gluing, but it pertains to an unique 
graph. Let Bi he a, D colored graph and select two vertices f 1 , wi € Bi. The contraction of B\ with the pair of 
vertices, denoted Bi/(wi,iJi) is the graph obtained from Bi by deleting the vertices vi and vi and reconnecting the 
lines touching them pairwise respecting the colorings. 

Again this operation can be performed in two steps. Consider the trace invariant associated to Bi and connect the 
two vertices vi and vi by a dashed line of color 0. Again this line will identify all the indices of the tensor associated 



to vi with the ones of the tensor associated to vi 
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SfU'if;^ , (3.18) 



as both vi and Vi are vertices in i3i. The contraction is represented graphicahy in figured] 
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Figure 4. Grapliical representation of the contraction of a graph. 

The contraction preserves the colorabihty. Note that the graph B/{vi,Vi) can potentially be disconnected (see 
figure S]). We will denote its connected components [B/(i'i,ui)] . Moreover, some of these connected components 
can consist in an unique line (as it is the case in figure H]). In this case they are not strictu sensu D colored graphs, 
but consist in exactly one line (with some color) which closes onto itself. This happens for every line which connects 
directly vi and ^;i in B. Such a line brings a factor N. At the level of the incidence matrix, one builds from e^g the 
incidence matrix 

v^v^.v^vr ef^/(''-^^) = ^ U (e% D e%) . (3.19) 

and the connected components with no vertices of course do not have an incidence matrix. 

These two operations encode the changing of an observable when adding a line of color 0. If the line of color in 
Q is a. tree line (joins Bi with some B2), then the observables Bi and B2 are glued. If on the other hand ?" is a loop 
line (starts and ends on the same Bi), then the observable is contracted. 

Geometrically these two operations have the following interpretation. Set for now D = 3. Then Bi and B2 represent 
surfaces. The surface associated to Bi is obtained by associating a triangle with edges colored 1, 2 and 3 to each 
vertex of Bi. This induces a coloring of the points (vertices) of the triangle by pairs of colors: 12 is the point common 
to the edges 1 and 2 and so on. A line in Si represents the unique gluing of the two triangles corresponding to its 
end vertices which respects all the colorings (those of the edges and those of the points, see [S^] for more details). 
The gluing Bi *(t,i.B2) ^2 comes to choosing a triangle (corresponding to vi) on Bi, a triangle (corresponding to V2) 
on B2 and gluing the two surfaces along the triangles. This is represented in figure [5] on the left, where we depicted 
the simplest case of the gluing of two planar surfaces. The contraction is essentially the same thing, just that this 
time the two triangles belong to the same surface, as represented in figure [S] on the right. 



Figure 5. Gluing and contraction of surfaces for tensors of rank D = 3. 

Note that the topology of the surfaces is changed under these moves (in the example of figure [5] on the right a planar 
surfaces becomes a genus one surface). This should come as no surprise as the same happens for matrix models: the 
contraction of a loop leads to two loops. The gluing is just a graphical encoding of surgery on the surfaces. The 
contraction has a more involved topological interpretation, and it can lead to an increase of the genus (if the two 
triangles contracted do not share anything), a splitting of the surface into connected components (if the two triangles 
contracted share at least two vertices but no edges) or no change at all (if the two triangles contracted share edges) . 

The most important feature of the gluing and the contraction is the following. Consider two observables Bi and B2 
joined by two lines of color 0, ('yi,i'2) and (^2,^1) with wi,wi G Bi and 'y2,'i'2 G S2. The resulting observable can be 
obtained in two ways, either by gluing along (wi,W2) and contracting with respect to (w2,^i) or the reverse. As the 
end result is unique, we have 

Ul,Wl G Bi , V2,W2 e B2^ [Bi *(„i,S2) 'B2]/(?;2,Wl) = [-62 *(t,2/Bi) Si]/(ui,U2) . (3.20) 



In D = 2 dimensions the observables are just bi colored cycles. The gluing of two cycles of lengths p and q always 
leads to a cycle of length p + q thus the gluing reduces to the addition of the lengths (hence it is associative). The 
colored gluing of graphs is associative (provided one tracks the vertices at which it is made) only if both expressions 
[;B*(„,gj) Bi] *(v',v2) ^2 and S*(„gj) [Si *(v'.v2) ^2] are defined, that is if v' € Bi. Note however that if v,v' £ B, 
while [S*(^ pj)] *(v',v2) ^2 is defined, S*(,„ gj^) [Bi '*^[v',v2) ^2] is not. The gluing is the appropriate generalization of this 
addition to the D-colored graphs representing the observables of tensor models. In the spirit of the matrix model 
nomenclature, one should call the trace invariants observables "bubble observables" , and the SDE's we derive in the 
next section "bubble equations". 



B. Schwinger Dyson equations and the algebra of constraints 



We now derive the SDE's of tensor models at all orders in N. We subsequently translate them into constraints 
satisfied by the partition function. The constraints form a Lie algebra, generalizing to all orders in N the D-aiy tree 
algebra identified in [53] . 

Consider a D-colored graph Bi. We chose a vertex vi ^ B (and mark it). For any Bi and -Di G Bi the following 
trivial identity holds 

Y.JidTdf]^^(5,^^,l[Tn^ n r,„ e e-^°-^«) = ^ 

J2 Trf5,/(„,,«,)(r,r))-iV^-i^t6^(^Tre.^„,,^,g,(r,r))=0. (3.21) 

We denote the \p\ connected components of Bi/{vi,vi) by [Si/(?;i, wi)] , with 1 < |p| < D. If one of these connected 
components consist in a single line the associated trace is N. Thus 

\p\ 

Tr^,/(„,,,,)(T,T) = n ^'Un..,..,] (^'^) ■ (3.22) 

p=i '- -' p 

The SDE's translates into a set of differential operators acting on Z indexed by the observable Bi and the marked 
vertex vi 
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(3.23) 



and by convention the partial derivative is — iV^ if [Si/(ui,wi)] is formed by an unique line. The natural domain 
of the differential operators Cm-^!g-^\ is the set of invariant functions depending on the coupling constants tg. The 
constraints form a Lie algebra 



Theorem 1. The commutator of two differential operators Cm^y-^-\ and C/ig-^y-^\ 
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Proof: The proof is a straightforward computation. We start from the commutator 
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The first line evaluates to 



d 



1 



d 
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(3.26) 



Consider the bubble Bi *{v,v2} ^2- When reducing with respect to (ui,wi) with vi in Bi it will disconnect into 



Bl*(v^V2) '^2 



several connected components < 

connected components [Si/(wi,wi)] . The special one is <. 

S^'P' , V r ,, 1 = y^ =R J hence the first line is 
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and exchanging the sums over v and vi it becomes 
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where we relabeled vi by v' and we added and subtracted the terms with i;' G S2. Recall that by equation p.20p . if 
v,vi e Z5i and v' ,V2 G B2 then 



^1 *(f ,02) ^2 
hence the first two lines in eq. (|3.25p yield 
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The third line in equation p.25p writes 
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We separate the terms with v' G Bi from the terms with w' G B to get 
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Adding equations (I3.30|) with (|3.31|) wc obtain, relabeling some dummy indices 

['C(Bi,Ci),-C(B2,«2)J = 
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This Lie algebra admits a closed Lie subalgebra. The leading order observables, the melons [3g, are indexed by 
colored rooted D-ary trees T- It is easy to check that gluing of the observables Bi *{vi,v2) ^2 reproduces the gluing of 
their associated trees Ti*v T2, as defined in [S^. The melonic observables are closed under this composition (as the 
gluing of trees leads to trees) , hence the algebra indexed by D-ary trees identified in [Sg is a Lie subalgebra of the 
full constraints algebra. 

IV. CONCLUSION 

We have derived in this paper the SDE's of tensor models for a generic complex tensor at all orders in 1/iV. They 
translate into a Lie algebra of constraints obeyed by the partition function. The algebra is indexed by colored graphs 
and generalizes to all orders in 1/N the algebra indexed by D-ary trees of the leading order observables. The algebra 
indexed by D-aiy trees closes a Lie subalgebra of the full constraints algebra. 

The study of this algebra of constraints, primarily of its central extensions and unitary representations, is a pre- 
requisite for the full classification of the continuum limits of tensor models. The study of its representations would 
benefit from identifying various Lie subalgebras, and studying the their representations first. We already posses a 
candidate, the leading order algebra indexed by D-ary trees. Other, simpler subalgebras can readily be identified: for 
instance the Virasoro algebra itself is a subalgebra of the full constraints algebra jSOf (in fact, following the results of 
|5(j | , one can identify several distinct copies of the Virasoro algebra as subalgebras of the full constraints algebra). The 
continuum SDE's should be understood in some appropriate double scaling limit and the continuum operators should 
be identified [isl . l20j . Other aspects of the emergent continuous geometry like its effective spectral and Hausdorff 
dimensions must also be analyzed. Analytic control of the continuum limit is a prerequisite in order to use the random 
tensor models to investigate conformal field theories, statistical models in random geometry and quantum gravity in 
arbitrary dimensions. 
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